variables = who('A*');

x=0;
y=0;

for i=1:size(variables)
    
    hist = eval([char(variables(i))]);
    
    s = hist(:,1);
    ns = hist(:,2);
    
    maxSize = 500;
    %maxSize = max(s);
    
    s = s(s<maxSize*.12);
    ns = ns(s<maxSize*.12);
    
    s = s(s>maxSize/100);
    ns = ns(s>maxSize/100);
    
    
    xx = z(epss(probabilities(i),pc),s,sigma);
    yy = ns.*s.^tau;
    
    x = cat(1,x,xx);
    y = cat(1,y,yy);
    
    plot(x,y,'.k')
end

